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Infrared, optical, and ultraviolet spectropolarimetry of cosmological sources is used to constrain 
the pure electromagnetic sector of a general Lorentz-violating standard-model extension. The coef- 
ficients for Lorentz violation are bounded to less than 3 x 10 -32 . 
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Lorentz violation is a promising candidate signal for 
Planck-scale physics [[[]. For instance, it could arise in 
string theory |2[] and is a basic feature of noncommu- 
tative field theories ||. In quantum field theory at at- 
tainable energies, small violations can be incorporated 
into the standard model to yield a general Lorentz- 
violating standard-model extension Q. Its lagrangian 
consists of all possible observer Lorentz scalars formed 
from standard-model fields while allowing for coupling 
coefficients with Lorentz indices. All renormalizable and 
gauge invariant terms relevant at low energies are explic- 
itly known. 

The standard-model extension predicts definite exper- 
imental signals. In the fermion sector of the theory, var- 
ious experiments have bounded coefficients for Lorentz 
violation. However, relatively little is known experimen- 
tally about the implications of the standard-model ex- 
tension for the properties of light. In particular, no 
bounds have been placed on the CPT-even coefficients 
for Lorentz violation in the photon sector. In this work, 
we study these terms and use spectropolarimetry of cos- 
mological sources to obtain stringent bounds on Lorentz 
violation comparable to the best current limits in the 
fermion sector. 

Extracted from the standard-model extension, the 
Lorentz-violating electrodynamics can be written in 
terms of the usual field strength Fp„ = d^Ay — d^A^, 
which incorporates the electric field E and the magnetic 
field B. The relevant lagrangian terms are Hpf 
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\[k F ) RXliV F KX F^. 



(1) 



The second term is CPT even and Lorentz violating. Its 
coefficient (kF) K \^v is dimensionless. It has the symme- 
tries of the Riemann tensor and zero double trace, so it 
contains 19 independent real components. 

The modified inhomogencous Maxwell equations ob- 
tained from Eq. (nl) are 



d a F^ + {k F )^ aPl d a F^ = 0. 



(2) 



These are supplemented with the usual homogeneous 
Maxwell equations d p F pv = 0. For a plane electro- 
magnetic wave with wave 4- vector p a = (p , p) , we have 
F^ v (x) = F llv {p)e~ %PaX . The homogeneous equations 



ensure the absence of single magnetic poles, p ■ B = 0, 
and provide the usual Faraday law, p°B = px E. This 
can be used to eliminate B in the modified inhomoge- 
neous equations, yielding the modified Ampere law 0] 

M jk E k = (_ 5 ik p 2 _ pJp k _ 2(k F y^ k P[3 p J )E k = 0. (3) 

The matrix AF fc is real and symmetric, which can be 
shown to imply lossless propagation. The modified 
Coulomb law p 3 M^ k E k — follows from this equation, 
in parallel with the usual case. 

For nontrivial solutions to Ampere's law, we must re- 
quire c\etA'P k = 0. This condition provides the dispersion 
relation between p° and p. Some calculation shows that, 
to leading order in the coefficients (kp) K \ fJ-u , the solutions 
to the dispersion relation take the form 



p\ = (l + p±a) \p\ 



(4) 
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where p = —\k a a and o 2 = \{k a p) 2 — p" 
(kpT^V^Vv and p» = p^/N- The solutions (|) de- 
scribe behavior similar to that of spatial dispersion in an 
optically inactive but anisotropic medium. In particu- 
lar, the vacuum is birefringent: light typically propagates 
as two different independent modes. Substitution of the 
two solutions (|J) into the Ampere law (||) determines the 
electric fields E± and hence the polarization of the two 
modes. For each mode, the group velocity v g ± = Vp- p± 
typically differs from the phase velocity v p ± = p^p/p 2 . 

At leading order in (fci^KA/ii/, the fields E± are or- 
thogonal and each is also perpendicular to v g ±. The 
orthogonality implies that E± span the set of all possi- 
ble polarizations, and so the unit vectors e± = E±/\E±\ 
form a basis for polarization at this order. The electric 
field can then be decomposed as E{x) = (E+i+e~' lp + t + 
E-i-e~ tp - t )e lp ' x . Since the phase velocities of the two 
modes differ, their relative phase changes as a wave prop- 
agates through free space. The resulting change in the 
polarization state of the wave is determined by the rela- 
tive phase change 



A<t> = (p° + -p°_)t w 2ttAv p L/\ « AitaL / \, 



(5) 



where Av p is the difference in phase velocities, A is the 
wavelength, and L is the distance the radiation has trav- 
eled. The distance dependence implies that, for sources 
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at very large distances, tiny differences in the phase ve- 
locities may be observable. 

In recent years, high-quality spectropolarimetry of dis- 
tant galaxies at infrared, optical, and ultraviolet frequen- 
cies has been performed [p|-^5[. These observations cor- 
respond to L/X greater than 10 31 . We can therefore an- 
ticipate that measurements of polarization parameters of 
order 1 should provide sensitivity of order 10 -31 or bet- 
ter to components of (kp) K \ fJ-l/ . The inverse dependence 
of A(f> on wavelength is a special feature exploited here 
to obtain a definite bound on (kp) K \^ v . 

In spectropolarimetry, the quantitative measurement 
of polarization is usually described via Stokes parame- 
ters, defined in a frame where the 3-axis coincides with 
the direction of propagation JTJ] . Introducing unit vec- 
tors ei and e-2 along the 1- and 2-axes and the corre- 
sponding electric field components as E\ and E2, the 
Stokes parameters can be taken as 

2Re.Ei*£:2, 2Im^i*£;2) 
= s (1, cos 2x cos 2tp, cos 2xsin2?/>, sin 2%), (6) 

where x an( i V 1 are the usual polarization angles. For 
convenience, we normalize throughout to s = 1. In the 
present context, nonzero coefficients (kp) K \^ u leading to 
a nonzero a imply a finite phase shift Acj) in the radi- 
ation from a cosmological source, which in turn affects 
the Stokes vector s(ip, \) determining its polarization. 
To set a bound on Lorentz violation, we must therefore 
first establish quantitatively the relationships between 
(&f)kA MI ,, c, A<f>, and s(ip,x)- 

We begin our analysis by expressing a directly using 
the 19 independent components of (kp) K \^ v . Since a 2 is 
a quadratic form in (kp) K x^, we can choose 19 indepen- 
dent components kA, A, B = 1, . . . , 19, of (kp)^^ and 
can write cr 2 = S^fe^fe^. The 19 x 19 matrix Sab is 
symmetric and depends on the direction of propagation. 
Some calculation shows that there exists a direction- 
independent similarity transformation such that Sab 
takes the form of a 19 x 19 matrix with only a 10 x 10 
nonzero block. Therefore, Sab has rank 10, and only ten 
linearly independent combinations k a , a — 1, . . . , 10, of 
(fcpOfcA/ii/ appear in a. An acceptable choice for these ten 
combinations is 

k a = ((kp) 0213 , (k F ) 0123 , 
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(M 0313 ). (7) 



It now follows that a 2 — S a bk a k b , where S a b is sym- 
metric and direction dependent. The other nine linearly 



independent combinations of (kp) K \^„ play no role in 
birefringence, and bounding them will be the subject of 
a separate investigation. For definiteness, the reference 
inertial frame in which the k a are specified by Eq. (0) is 
chosen to be compatible with celestial equatorial coordi- 
nates, with the 3-axis aligned along the celestial north 
pole at equinox 2000.0 at a declination of 90°. The 1- 
and 2-axis are at a declination of 0° and a right ascen- 
sion of 0° and 90°, respectively. The goal is to bound the 
ten quantities k a defined in this frame. 

The form of S a b is cumbersome and is omitted here. A 
more convenient expression for a 2 can be obtained by cal- 
culating within a special inertial frame. The idea is to use 
observer rotation invariance to choose a 'primed' frame 
in which has leading-order form = (1; 0,0,1), to 
solve for the Lorentz scalar a 2 in this frame in terms of 
(^-fOkV' an< ^ then to use the rotation between the ce- 
lestial equatorial frame and the primed frame to express 
(kp)' KXjlv in terms of (kp) K Xfj,u- To match standard po- 
larimetric conventions, we choose the primed-frame basis 
vector e 3 to point from the source towards the Earth, 
while e[ points south [T^j . 

Solving the Ampere law (|^) in the primed frame gives 
p = i(fe m + fc' 22 ) and a 2 = (fc' 12 ) 2 + ±(fc' n - fc' 22 ) 2 . 
This form for a 2 suggests defining an angle £ such that 
fc' 12 = crsin£ and §(fc' n - A:' 22 ) = <jcos£. The angle £ 
determines the minimal linear combinations of (kp)^^ 
relevant for polarimetry of a specific source. Note that £ 
is frame dependent, unlike p and a. 

At leading order, the polarization basis vectors in 
this frame obey e± oc (sin£, ±1 — cos£, 0) and are lin- 
early polarized. The corresponding Stokes vectors are 
s± = ±(cos£,sin£,0). The propagation from the source 
to the Earth generates a relative phase change Acj) spec- 
ified by Eq. (g). The corresponding effect on the Stokes 
vector s(if>, x) f° r the radiation can be regarded as a ro- 
tation by A(f> about the basis vector s+. This typically 
changes both if) and \. 

The change in polarization depends not only on the co- 
efficients (kp) K \ tll/ , the wavelength A, and the distance to 
the source L, but also on the initial polarization. For cos- 
mological sources, there is usually no way to determine 
independently the polarization produced at the source. 
We adopt instead the reasonable assumption that the 
source polarization is constant over the relatively narrow 
band of wavelengths being considered for each source. 
The quantity of interest is then the change in relative 
phase 8(f) — Ana (L/X — L/Xq) between a wavelength A 
and a reference wavelength Ao- 

The rotation of the Stokes vector can be expressed via 
a Mueller matrix m? k , often used to describe the effects 
of filters and polarizers on light |l6| . The change in the 
Stokes vector is given by s^ip.x) — m (<^>) s (VtojXo)) 
where "0o, Xo are reference polarization angles. Some al- 
gebra reveals that 
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cos 2 £ + sin 2 £ cos < 
(50)= I sin £ cos £(1 — cos S 
— sin £ sin 56 



sin£ cos£(l — cos 6 
sin 2 £ + cos 2 £ cos < 
cos £ sin 50 



sin£ sin Si 
- cos £ sin < 
cos 86 



(8) 



The polarization angle x, which measures the amount 
of circular polarization, is unavailable in most published 
literature on astronomical spectropolarimetry. We there- 
fore focus attention here on the change Sip = ip — ipo in 
the polarization angle ip. After some algebra, we find 



Figure 1 shows the change in the minimized distribu- 
tions from their least value for the source 3CR 68.1. The 



Sip = h tan 



_ x sin £ cos Co + cos £ sin Co cos{86 ~ 6 ) 
cos £ cos Co — sin £ sin Co cos(S6 — 6q) 



, (9) 



tan 1 (tan2xo/sin£), 



where £ = £ — 2ipo and 6a 
Co = cos~ 1 (cos2xocos£). 

To obtain a bound on the coefficients (/cf^a^i/, our 
strategy is to fit Eq. (^|) to polarization measurements 
of distant sources at multiple wavelengths. Since a sin- 
gle source only constrains a and £ in the particular 
primed frame appropriate for that one source, a mea- 
surement involves only a two-dimensional surface in the 
ten-dimensional coefficient space. It follows that at least 
five different sources are required to make a complete 
measurement of all the coefficients k a . 



Source 



IC 5063 
3A 0557-383 
IRAS 18325-5925 
IRAS 19580-1818 
3C 324 
3C 256 
3C 356 

F J084044.5+363328 

F J155633.8+351758 

3CR 68.1 

QSO J2359-1241 

3C 234 

4C 40.36 

4C 48.48 

IAU 0211-122 

IAU 0828+193 



L (Gpc) 



0.04 
0.12 
0.07 
0.13 
1.69 
1.92 
1.62 
1.71 
1.82 
1.70 
1.48 
0.55 
2.02 
2.04 
2.04 
2.08 



10 M L/\ 



0.56 - 2.8 
2.2 - 8.4 
1.0 - 4.9 
1.8 - 9.1 

58 - 130 
70 - 140 
57 - 120 
62 - 120 
67 - 110 

59 - 130 
87 - 90 
51 - 75 

73 - 160 
75 - 160 

74 - 160 
78 - 160 



log 



10 ' 



-30.8 
-31.2 
-31.0 
-31.0 
-32.2 
-32.2 
-32.2 
-32.2 
-32.2 
-32.2 
-31.1 
-31.7 
-32.2 
-32.2 
-32.2 
-32.2 



Rcf. 



[S| 
[71 
[71 
[7] 
[3[ 
[9| 
[10| 

iyi 
iyi 

[12] 
[13] 

\W\ 
[15| 
[15| 
[15| 

Pi 



Tabic 1. Source Data. 

Since the sensitivity to k a is inversely related to wave- 
length, polarimetry at shorter wavelengths yields better 
bounds. To maximize the sensitivity, we restrict atten- 
tion to a chosen sample of 16 sources with published po- 
larimetric data well suited to our purpose. Table 1 lists 
this sample, which spans wavelengths in the range 400- 
2200 nm. For each source, we choose ipo as the mean 
polarization angle and use Eq. (JTJ) to create a \ 2 distri- 
bution for Sip as a function of a, £, Ao, Xo- For each a, £, 
we minimize with respect to Ao and xo, which produces 
values consistent with expectations. 
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50 is displayed in 



the £-log a plane. This corresponds to confidence level 
(100 — 10 _9 )%, which we regard as sufficient for a defi- 
nite bound. Only —90° < £ < 90° need be considered, 
since Sip is symmetric under £ — * £ ± 180°. The gen- 
eral shape of this plot is common to all sources. Each 
source eliminates some region of coefficient space away 
from a = 0, £ = 0,±90°, which are the only configu- 
rations for which the theory predicts no change in the 
polarization angle. In principle, by combining data from 
multiple sources it would be possible to eliminate the re- 
gions near £ = 0, ±90° for all the sources. However, in 
the present context it suffices to make the reasonable as- 
sumption that the true signal does not lie in these regions 
for all the sources. Under this assumption, we have cho- 
sen conservative bounds on a for each source, listed in 
the fourth column of Table 1. For the source 3CR 68.1, 
the bound is shown as a horizontal line across Fig. 1. 
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FIG. 1. Contour with A* 2 = 50 for the source 3CR 68.1. 

To estimate the constraint on k a that results from com- 
bining these bounds, we assume that the data for each 
source are consistent with a = 0. The bounds on a can 
then be regarded as conservative estimates of the error 
So in a null measurement of a. We then create a second 
X 2 distribution, x 2 — J2j ( a j) 2 /(^ a j) 2 > where j ranges 
over the 16 sources. Constant values of this distribution 
correspond to ten-dimensional ellipsoids in the coefficient 
space. Taking the magnitude |fc a | = Vk a k a of k a as the 
variable of interest and minimizing \ 2 with respect to the 
other nine degrees of freedom produces a constraint on 
|fc a |. Our conservative estimates yield |fc a | < 3 x 10~ 32 
at the 90% confidence level. 

As a check on this constraint, we have performed a 
grid search on the ten-dimensional space of coefficients 
k a . Fixing the magnitude of the ten-dimensional vector 
k a in Eq. (Q) leaves a nine-dimensional sphere. For a 
given point on the sphere, a and £ for each source can 
be found and the value of x 2 obtained using the mini- 
mized \ 2 plots. Summing over sources yields a net x 2 
minimized with respect to all parameters other than the 
coefficients k a . A grid search on the nine-sphere with 
points separated by about 5 degrees of arc was used to 
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find the minimum net \ 2 ■ Repeating for different values 
of \k a \ gives results confirming our bound. To check our 
procedure directly, we also generated synthetic polariza- 
tion data with a known faux Lorentz- violating signal and 
verified that our procedure correctly extracts it. Note 
that the absence of a signal emerging from our analysis 
indicates that systematic effects are irrelevant. To gener- 
ate a false signal, these would need to mimic the expected 
direction and wavelength dependence. 

No previous bounds on the coefficients (kp) K \ tll/ ex- 
ist. However, constraints on different Lorentz-violating 
coefficients in the fermion sector of the standard-model 
extension have been obtained from studies of neutral- 
meson oscillations |l^-|2l|l , comparative tests in Penning 
tra ps |2^ -|25|l , spectroscopy of hydrogen and antihydrogen 
p(i| . [27| , measurements of muon properties ]2S|jffi| , clock- 



comparison experiments 30 - 32 1 , observations of the be- 
havior of a spin-polarized torsion pendulum (3^J|^], and 
studies of the baryon asymmetry |35f| . The constraint re- 
ported here is comparable to the best of these existing 
limits, presently a few parts in 10 . 

An improved bound could be obtained as more high- 
quality data become available, particularly if measure- 
ments of the polarization angle x could be incorporated 
in the analysis. Also, the sensitivity to inverse wave- 
length implies that spectropolarimetry of cosmological 
sources at frequencies above the ultraviolet band would 
be of interest. The technology to undertake X-ray po- 
larimetry of cosmological sources is presently being de- 
veloped p6| , which suggests an improvement of several 
orders of magnitude may eventually be attainable. 

This work was supported in part by DOE grant DE- 
FG02-91ER40661 and by NASA grant NAG8-1770. 
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